# -*-coding:utf-8-*-

# 执行pssm文件

import os



for root, dirs, files in os.walk(r"fasta2data"):
    for file in files:
        print(root)
        print(os.path.join(root, file))
        os.system(
            rf'psiblast.exe -query .\fasta2data\{file} -db D:\install\NCBI\blast-2.12.0+\db\sirt\sirt -comp_based_stats 1 -inclusion_ethresh 0.001 -num_iterations 3 -out_ascii_pssm .\taipssm\{file[:-6]}.pssm')
